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Abstract 

The Tracy- Widom distribution functions involve integrals of a Painleve II function starting from 
positive infinity. In this paper, we express the Tracy- Widom distribution functions in terms of integrals 
starting from minus infinity. There are two consequences of these new representations. The first is the 
evaluation of the total integral of the Hastings-McLeod solution of the Painleve II equation. The second 
is the evaluation of the constant term of the asymptotic expansions of the Tracy- Widom distribution 
functions as the distribution parameter approaches minus infinity. For the GUE Tracy- Widom distri- 
bution function, this gives an alternative proof of the recent work of Deift, Its, and Krasovsky. The 
constant terms for the GOE and GSE Tracy- Widom distribution functions are new. 

1 Introduction 

Let Fi (x) , F2 (x) , and F^(x) denote the GOE, GUE, and GSE Tracy- Widom distribution functions, respec- 
tively. They are defined as [221 US] 

F 1 (x) = F(x)E(x), F 2 (x)=F(x) 2 , F 4 (x) = \\e{x) + i W 

where 

F(x) = exp(-^ R(s)dsj, = exp - i ^ q(s)dsj. (2) 

Here the (real) function q(x) is the solution to the Painleve II equation 

q" = 2q 3 + xq, (3) 

that satisfies the boundary condition 

q(x) ~ M(x), x — > +00. (4) 
Recall [1] that the Airy function Ai(x) satisfies Ai"(x) = xAi(x) and 

A •/ \ 1 -2 3/2 

MX) ~ 2^/4 e ' ^+°°- ( 5 ) 
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There is a unique global solution q(x) to the equation with the condition (fj| (the Hastings-McLeod 
solution) [15]. The function R(s) is defined as 

/>oo 

R(x) = / (g( S )) 2 ds. (6) 



By taking derivatives and using and ((D) (see. for example, (1.15) of [55] and (2.6) of [3]), the function 
R(x) can also be written as 

R(x) = (q'(x)) 2 -x(q(x)) 2 -(q(x))\ (7) 
Integrating by parts, F(x) can be written as 

F(x) = exp ( - J / ( S -.T)(g( S )) 2 d S V (8) 



which is commonly used in the literature. 

Notice that ([2|) involves integrals from x to positive infinity. The main results of this paper are the 
following representations of F(x) and E(x), which involve integrals from minus infinity to x. 

Theorem 1.1. For x < 0, 

**> = 2"«e«'<-^exp{l £ (hM - ^ + i-)*}, (9) 
where £(z) is i/ie Riemann-zeta function, and 

^) = ^7re-^ |x|3/2 exp{i f - JM'Ul. (i ) 



21/4^ ^12_ 

Remark 1. TTie formula @ afeo follows from the recent work JSjj of Deift, Rs, and Krasovsky. See subsec- 
tion below for further discussion. 

The integrals in © and (JTUJ) converge. Indeed, it is known that [TH [T2] 

1 73 . 10219 . ,_ 12 , 



= V - I 1 + &? - 128^ + 1024^ + ° (N V ' ^ ^ (U) 

This asymptotic behavior of q was obtained using the integrable structure of the Painleve II equation (see, for 
example, |16j). The coefficients of the higher terms in the above asymptotic expansion can also be computed 
recursively (see for example, Theorem 1.28 of [12]). For R(x), |Zj) and (fTT|) imply that 

, v x 2 ( 1 9 128 , 19 A , . 

We now discuss two consequences of Theorem 11.11 

1.1 Total integrals of q(x) and R(x) 

Comparing with (|2|), Theorem 1 1.1 1 is equivalent to the following. 

Corollary 1.2. For c < 0, 



r°° r c / i 1 \ 1 11 

J R(y)dy + J_ f J R(j / )-- y 2 + _j^ = __ log2 _ C / ( _ 1) + _| c |3 + _ log | c | 



(13) 



1 , n , V^, ,3/2 



q(y)dy+ I \q{y)~\l^)dy = -\o g 2+^\c\ i /'. (14) 
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These formulas should be compared with the evaluation of the total integral of the Airy function 

ID: 

/OO 
Ai{y)dy = 1. (15) 
-oo 

Recall that the Airy differential equation is the small amplitude limit of the Painleve II equation. Unlike 
the Airy function, R(x) and q{x) do not decay as x — > — oo, and hence we need to subtract out the diverging 
terms in order to make the integrals finite. 



1.2 Asymptotics of Tracy- Widom distribution functions as x — ► — oo 

Using formulas ([2]) and (fT2|) . Tracy and Widom computed that (see Section l.D of [24]) as x — > — oo, 

F 2 (x) = T2^^(l + ^ + 0(M- 6 )) (16) 
for some undetermined constant t%. The constant r 2 was conjectured in the same paper [24] to be 

T2 = 2 l/24 e C'(-l)_ (17) 

This conjecture (|17|1 was recently proved by Deift, Its, and Krasovsky [8j. In this paper, we present an 
alternative proof of (fTT|) . Moreover, we also compute the similar constants t\ and T4 for the GOE and GSE 
Tracy- Widom distribution functions. The asymptotics similar to (|16[) follow from {!]) and f| 1 lj) : as x — > —00, 



Fifi) = Ti j — ttttt 1 F r + 0(b| _d ) , (18) 

U ; M 1 / 16 V 24V2|x| 3 /2 Vl 1 7 ' 

F4(x) = T4 — ™ — ( 1+ w^ +0(N }) 



|1/16 

Using (jTTJ) and Q12[) . Theorem 11.11 implies the following. 
Corollary 1.3. ^4s a; — » 00, 



and 



Hence 



Fix) = 2 V« e K(" 1 )^ (l + 3 + 0(|,|- 6 ) ) 



= ^Tre"^ 1 "' 372 ( 1 ^ + 0(|x|- 3 ) ) . (21) 

V ^ 2V4 ^ 24V2>| 3 / 2 Vl ' 7 1 ; 

n =2- 11 / 48 e^'(- 1 ), T 2 = 2 1 / 24 e^- 1 \ r 4 = 2- 35 / 48 e K'(-D. (22) 



Conversely, using (fTTj) and (fl^j) . this Corollary together with (f5|) implies Corollary 11.21 and hence 
Theorem 11.11 

This is one example of so-called constant problems in random matrix theory. One can ask the same 
question of evaluating the constant term in the asymptotic expansion in other distribution functions such as 
the limiting gap distribution in the bulk or in the hard edge. For the gap probability distribution in the bulk 
scaling limit which is given by the Fredholm determinant of the sine-kernel, Dyson [13] first conjectured the 
constant term for (3 = 2 in terms of C'(~ 1) using a formula in an earlier work [57] of Widom. This conjecture 
was proved by Ehrhardt [14] and Krasovsky [20] , independently and simultaneously. A third proof was given 
in [9]. The constant problem for f3 = 1 and (3 = 4 ensembles in the bulk scaling limit was recently obtained 
by Ehrhardt [15]. For the hard edge of the /3-Laguerre ensemble associated with the weight x m e~ x , the 
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constant was obtained by Forrester [TH] (equation (2.26a)) when m is a non- negative integer and 2//3 is a 
positive integer. 

The above limiting distribution functions in random matrix theory are expressed in terms of a 
Fredholm determinant or an integral involving a Painleve function. For example, the proof of [5] used the 
Frcldhom determinant formula of the GUE Tracy- Widom distribution: 

F 2 (x) = det(l-A SB ), (23) 

where A^ is the operator on L 2 ((x, oo)) whose kernel is 

u — V 

In terms of the Fredholm determinant formula, the difficulty comes from the fact that even if we know all the 
eigenvalues Xj(x) of A^, we still need to evaluate the product n^. 1 (l — Xj(x)). When one uses the Painleve 
function, one faces a similar difficulty of evaluating the total integral of the Painleve function. 

We remark that the asymptotics as x — > +00 of F(x) and E{x) (and hence Fp(x)) are, using Q, 

F(x) = 1 
E{x) = 1 



1 



35 



327TX 3 / 2 V 24z 3 / 2 

_ 2 _3/2 



+ 0(x~ i ) 



1 



41 

48x 3 / 2 



+ 0(x" J ) 



(25) 
(26) 



1.3 Outline of the proof 

The Tracy- Widom distribution functions are the limits of a variety of objects such as the largest eigenvalue 
of certain ensembles of random matrices, the length of the longest increasing subsequence of a random 
permutation, the last passage time of a certain last passage percolation model, and the height of a certain 
random growth model (see, for example, the survey [2 1 j ) . Dyson |13j exploited this notion of universality 
to solve the constant problem for the sine- kernel determinant. Namely, among the many different quantities 
whose limit is the sine-kernel determinant, he chose one for which the associated constant term is explicitly 
computable (specifically, a certain Toeplitz determinant on an arc for which the constant term had been 
obtained by Widom |27j ) , and then took the appropriate limit while checking the limit of the constant term. 
However, the rigorous proof of this idea was only obtained in the subsequent work of Ehrhardt [14j and 
Krasovsky |20j . In order to apply this idea for F%{x), the key step is to choose the appropriate approximate 
ensemble. In the work of Deift, Its, and Krasovsky [8], the authors started with the Lagucrre unitary 
ensemble and took the appropriate limit while controlling the error terms. In this paper, we use the fact 
that Fp(x) is a (double-scaling) limit of a Toeplitz/Hankcl determinant. 

Let D n (t) denote the n x n Toeplitz determinant with symbol f{e ie ) — e 2tcos W on the unit circle: 



D n (t) = dot (J- J e 2tcos V (j - fe)e d(9 



0<j,k<n-l 



1 



(27r) n n! 



(27) 

3 2t£? =1 co 8 e, \e ie " -e^fUdOj. 

l<k<£<n j=l 



Note that some references (e.g. [7]) define D n as an (n + 1) x (n + 1) determinant, whereas others (e.g. 
[5]) use our convention. In studying the asymptotics of the length of the longest increasing subsequence in 
random permutations, in [2], the authors proved that when 

n= [2t + xt 1/s ], (28) 

as t — > 00, 

e-*D n (t)->F 3 (x). (29) 
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The idea of the proof of (|2U)) in [3J is as follows. The Toeplitz determinants are intimately related to 
orthogonal polynomials on the unit circle. Let Pj(z) = KjZ J + • • ■ be the orthonormal polynomial of degree 



j with respect to the weight e 2 * cos9 |i : 



p j (e»)p k {e*)e tt ' sa "'— = 5 jk for j, k > 0. (30) 

If Kj > then pj is unique. We denote by tTj(z; t) = TTj(z) the monic orthogonal polynomial: Pk(z) = Kk^j(z). 
Then (see, for example, |23| ) the leading coefficient Kj = Kj(t) is given by 



As the strong Szego limit theorem implies that D n (t) — > e* 2 as n — * oo for fixed t, the left-hand-side of (|29p 
can be written as 

e~ t2 D n (t) = f[^f ) =f[<(t)- (32) 

q—n ^"^ ^ ' q—n 

The basic result of [2] is that 

^W-i-ffTT' q=[2t + yt 1 / 3 } (33) 

for y in a compact subset of R. (In [2J, the notations v(x) = —R(x) and u(a;) = —q(x) are used.) Hence 
formally, as f — > oo with n~2t + xt 1 / 3 , 

log (e-* 2 £>„(*)) = f] log («»(t)) ~ t 1 / 3 log (l - ~ - f° R(y)dy = logF 2 (x). (34) 

q=n x » ' 31 

The first step of this paper is to write, instead of ([32]) . 

Here Do(t) := 1. Then formally, we expect that as f — » oo with n = 2t + xt 1 ^ 3 , (|35|) converges to an integral 
from — oo to x. For this to work, we need the asymptotics of K q (t) for the whole range of q and t such that 
1 < q < 2t + xt 1 / 3 as t -> oo. 

It turns out it is more convenient to write, for an arbitrary fixed L, 

e"* 2 D n (t) = e"* 2 D L (t) f[ -M> - = D L (t) f[ (36) 

q =L + l F i + lVlW 

We introduce another fixed large number M > and write 

[H-Mt 1 ' 3 -!} [2t+xt 1/3 ] 

log(e- t2 D n ) = -t 2 + log(D L )+ l0 S(«7-x)+ E !°g(«7-i)- ( 37 ) 

, ^ L+1 , ?=[2t-M*V3] 

v v ' * „ 

Airy part Painlevc part 

Since L and M are arbitrary, we can compute the desired limit by computing 

lim lim log(e"* 2 £)„(i)), n = [2t + xt 1/3 }. (38) 

L,M— too t^oo 
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From (|33p . the Painlcvc part converges to a finite integral of R(y) from y = —M to y = x as t — > oo. 
For the Airy part, we need the asymptotics of /c g (i) for L + 1 < g < 2t — Mt 1 ' 3 as t — > oo for fixed 
L, M > 0. The paper [5] obtains a weak one-sided bound of K q (t) for et < q < 2t — Mt 1 ' 3 as i — > oo, where 
e > is small but fixed. The technical part of this paper is to compute the leading asymptotics of K q (t) in 
L + 1 < q < 2t — Mt 1 / 3 with proper control of the errors so that the Airy part converges. The advantage 
of introducing L is that we do not need small values of g, which simplifies the analysis. The calculation is 
carried out in Section [3] Finally, for the exact part, the asymptotics of Di,(t) as t — > oo are straightforward 
using a steepest-descent method since the size of the determinant is fixed and only the weight varies. The 
limit is given in terms of the Selberg integral for the L x L Gaussian unitary ensemble, which is given by a 
product of Gamma functions, the Barnes G-function. The asymptotics of the Barnes G-function as L — ► oo 
are related to the term C'( — 1) ( sce P%)) below). The computation is carried out in Section^ 

Now we outline the proof of the formula (fTU|) for E(x). In the study of symmetrized random 
permutations it was proven in [H [5] that, in a similar double scaling limit, certain other determinants 
converge to -Fi(x) and F±(x). But it was observed in [?J[5] that these determinants can be expressed in terms 
of K q (t) and iTq(0;t) for the same orthonormal polynomials (|30[) above. Hence by using the same idea for 
D n (t), we only need to keep track of 7r g (0;i) in the asymptotic analysis of the orthogonal polynomials. See 
Section [5] below for more details. 

This paper is organized as follows. In Section [21 the asymptotics of the exact part of (|37| arc 
computed. We compute the asymptotics of the Airy part in Section [3] The proof of for F(x) in 
Theorem 1 1.1 1 is then given in Section 0] The proof of (|10[) for E(x) in Theorem II. II is given in Section [5l 

While we were writing up this paper, Alexander Its told us that there is another way to compute 
the constant term for E(x) using a formula in [5]. This idea will be explored in a later publication together 
with Its to compute the total integrals of other Painleve solutions, such as the Ablowitz-Segur solution. 

Acknowledgments. The authors would like to thank P. Deift and A. Its for useful communications. The 
work of the first author was supported in part by NSF Grant # DMS-0457335 and the Sloan Fellowship. The 
second and third authors were partially supported by NSF Focused Research Group grant # DMS-0354373. 



2 The exact part 

We compute the exact part of ([37]) . From equation (|27|) . 

a 2t Ef =1 cos g j TT I JOk _ je e | 2 

(2w) L L\ J<_-„- 



D ^) = 7^Tn[ e*£k~* I] \e^-e^\ 2 Y[d9 r (39) 



i<k<e<L j=i 

Following the standard stationary phase method of restricting each integral to a small interval —e<0<e 
and expanding e and e 2tcos9 in Taylor series, D^t) is approximately 

v ; •'L-e.eJ l<k<l<L j=l 

as t — * oo. By extending the range of integration to R L , we obtain 

}^ D ^i(2^ D ^) =i « ( 41 ) 

where 

DT m (t) = 1 / e-^U e > J] - II dB r ( 42 ) 

^- J[-oo,oo]i l<fe<^<„ j=l 
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This integral is known as a Selberg integral and is computed explicitly as (see for example, |22j . equation 
(17.6.7)) 



^Hcrm 



(*) = 



r L/2 



L-l 



2 L(L-l)/2 t L*/l 



11*' = 



L/2 



q=0 



2 L(L-l)/2 t L2/2 



G(L + 1), 



(43) 



where G(z) denotes the Barnes G-function, or double gamma function. Some properties of the Barnes 
G-functions are (see, for example, [2"6l IS] ) 



G(z + l)=T(z)G(z) 
G(l) = G(2) = G(3) = 1. 

l0gG (l) =^l°g2- Jlog^C'(-l), logGQ) =-Llog2+ilog7r+^C'(-l) 

G{n) = l!2!---(n-2)!, n = 2, 3, 4, • • • . 

z 2 3 z 1 

logG(z + 1) = — logz - -z 2 + - log(27r) - — logz + C'(-l) + O I 



Therefore, 



^ Urn \ \og{D L ) - j 2Lt-^-log(2t) + _ L\ iogL - l L 2 + C '(-l) > 



(44) 
(45) 
(46) 
(47) 
(48) 

(49) 



3 The Airy part 

The main result of this section is Lemma 13.101 which computes 



lim 

L,M— >oo t— kx> 



the Airy part of ([57)) . We use the notation 



[2t-Mt 1/3 -l] 

iog(^-i), 



_ 2i 
/ 9 ' 



(50) 



(51) 



It is well known that the leading coefficients n q 2 of orthonormal polynomials can be expressed in 
terms of the solution of a matrix Riemann-Hilbert Problem (RHP) (T7J. We start with the RHP for 
defined in Section 6 (p. 1156) of [2], which is obtained through a series of explicit transformations of the 
original RHP for orthogonal polynomials. For notational ease, we drop the tildes Let 6 C be defined such 
that < 6 C < 7r and sin 2 = • ^ or 1 ln ^ ne regime L + 1 < q < 2t - Mt 1 / 3 - 1, we have 7 > 1. Define 
the contours C\ = {e 10 : 6 C < \6\ < ir} and G2 = {e 10 : < \8\ < 8 C } with the orientations given as in 
Figure 1(a) Also define the contours Gi n and G ou t as in Figure 1(a) Let S^ 5 ' = G\ U G2 U G; n U G ou t. Now 
let m^(z) = mS 5 ' {z;t,q) be the solution to the following RHP: 



to' 5 ' is analytic in z € C\Es 



•?>(*) 
(5) = / as z 



,(5) 



z)t/ 5 )(.z) for z G £5 



(52) 



00 
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where the jump matrix = v^(z; t, q) is given by 



v {5) (z) = i 



Here 




where £ = e c and the branch is chosen to be analytic in C\C2 and v(s — £)(s — £) 



(53) 



(54) 

-s for s — » oo. Then 




(a) The contour S< 5 > for m( 5 >. 



(b) Introduction of and 



(c) The contour S for m^' . 



Figure 1: Contours used in the definition of the RHP for and rrS R \ 



(see (6.40) of g]) 



and 



- _ P 9(-7+log7+l) Trl ( 5 ) 



m^(0) 



(55) 



7r,(0) = (-l)Xr(0). (56) 

We analyze the solution to' 5 ' to this RHP for the regime L + l<q<2t — Alt 1 / 3 — 1 as t — > oo. Our 
analysis builds on the work of [2] and makes two main technical improvements. The first is that the paper 
[2] only considered the regime when et < q. Hence q necessarily grows to infinity. In this work, we allow q 
to be finite. The second is that we compute a higher order correction explicitly to the asymptotics obtained 
in [2] . This higher-order correction contributes to the sum (|50|) . In [2] , only a one-sided bound of a similar 
sum was obtained. We merely outline the analysis for the parts that overlap with the analysis of [2J. 

From the construction of a in g] we have that le - "' 2 -'! < 1 for z <G C\ and | ? | < 1 for z £ 
C; n U C ou t- If we formally take the limit of our jump matrix v^ 5 ' as q — ► oo the jumps on the contours C m 
and C ou t approach the identity matrix and the jumps on C\ and Ci approach constant jumps. This limiting 
RHP is solved explicitly by 



,(5,00) 



(z) 



kiP+P' 1 



(57) 



/ _ f \l/4 

where /3(z) — ( z f^-i J , which is analytic for z G C \ C2 and /3 — > 1 as z — > 00. Note that 



(6,< 
'21 



,(5/ 



'(0) 



2/ 



2/ 



(58) 



However, the convergence of the jump matrix z/ 5 )(z) is not uniform near the points £ and £, since 
a(£) = a(£) = 0. Therefore a parametrix is introduced around these points. For fixed S < y^j, define 

Ot = {z:\z-Z\<5\Z-l\}, O i ={z:\z-t\<8\£-Z\}. (59) 

Note that the diameter of is of order and varies as t and q vary. The diameter approaches as 

7 — ► 1 or 7 — > 00, which happens when q is close to 2t — Mt 1 ^ 3 or L + 1, respectively. However, the point is 
that in the regime L + l<q<2t— Alt 1 / 3 — 1, the diameter of cannot shrink "too fast." Therefore, the 
usual Airy parametrix for a domain of fixed size still yields a good parametrix for the RHP in the regime 
under consideration. The case when 7 — > 1 "slowly" was analyzed in [2] for the leading asymptotics of mS 5 ' . 
In this section, we also analyze the case when 7 — > 00 "slowly," and also improve the work in [2] to obtain a 
higher-order correction term. 

Orient the boundary of both and in the counterclockwise direction. Now as in [5] (see also 
|12| ) for z G 0£ \ £( 5 ) define the matrix- valued function m p as 



m p (z) = 
where lo = e 27 ™/ 3 and 



1 



lWwf?«(*)Y(^i 




qa(z)o- a 



(60) 



= < 



Ai(s) Ai(u> 2 s) 

Ai'(s) w 2 Ai'(w 2 s) y 
Ai(s) Ai(u> 2 s) 

Ai'(s) w 2 Ai'(w 2 s) y 

Ai(s) -w 2 Ai(ws) 

Ai'(s) -Ai'(ws) 

Ai(s) -w 2 Ai(cjs) ' 

Ai'(s) -Ai'(ws) 



< arg(s) < f , 



e e 



e s 




2 TV 
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< arg(s) < 7r, 
7T < arg(s) < 
^ < arg(s) < 2tt. 



(61) 



We can define m p (z) = m p (2) for z€^by {60]). For z ^ T,^ U (£\ U let m p (z) = m^^\z). It is 
shown in [2] that m p then solves a RHP that has the same jump conditions as to' 5 ' on the contour C2 as 
well as on T,^ n O, where we define O = U Oj-. 

Define i?(z) = m^{z)m- x {z). Then i? solves a RHP on S = dO U ((Ci U C in U C out ) n O c ) with 
jump = m p _u( 5 ^i;~ 1 TOp_ 1 :. Explicitly, the jump matrix Vr is given by 



vr = I + < 




ze (s( 5 )nO)uC 2 , 



(62) 



z g (Cin u Cout) n o c , 

z G dO, 
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where v q p is given explicitly in Lemma 13-11 below . and the matrix v% is defined as tif 



vr 



z e dO. Since to^O) = R(0)m p (0) = i?(0)m( 5 '°°)(0), we have m^O) 



I —vl? for 

qa H 



^(0) = -1(^(0) - v^TiiuCO)). Therefore 



^(i? 22 (0)-V^ 



%/7 



l<?(-7+log 7+1)^(0) _ v /^3li? 21 (0)) 



li?2i(0)) and 



(63) 



and 



7T 9 (0) 



-1)? 



(i?i 2 (0)- Vt^T^hW). 



(64) 



1 qa 
qa tt 



In [2J, for z G <90, the jump matrix vr is approximated by the identity matrix / and the terms 
Vft are treated as an error (for the case when et < q). For our purpose, we need to compute the 
contribution from the next order term -^Vr* explicitly. 

It will be shown in the following subsections that for any e > 0, there are Lq and Mq such that 
for fixed L > Lq and M > Mq, there is tg — to(L,M) such that \\vr — I\\l°°C£) < e f° r au * — an d 
L + l<q<2t- Alt 1 / 3 - 1. This was shown in [2] for e x t < q < 2t - Alt 1 / 3 - 1. Then we proceed via 
the standard Riemann-Hilbert analysis as, for example, in [2j. Let €(f)(z) := -Jh f„ ^zr^ds be the Cauchy 
operator defined for z ^ S. For z e E, is defined as the nontangential limit of £(f)(z') as z' 

approaches z from the right-hand side of S. Define the operator £r(/) = <£-(f(vR — I)) for / e £ 2 (£) and 
the function /i = / + (1 — <£r)~ 1 <ZrI . A simple scaling argument shows that Cr is a uniformly bounded 
operator for L + 1 < q < 2t - A/t 1 / 3 - 1. Since the supremum norm of vr — I can be made as small as 
necessary, we find that for L and M fixed but chosen large enough, (1 — Cr) -1 is a bounded L 2 operator 
with norm uniformly bounded for t sufficiently large for all q such that L + 1 < q <2t — Alt 1 / 3 — 1. By the 
theory of Riemann-Hilbert problems, 



R{z) - I 



1 

2ttI 



I^{S)(VR - I) 



ds. 



(65) 



Define the contours as the part of £ in the upper-half and lower-half planes, respectively. That 

E ± =Sn(±3(z)>0). (66) 



Also define 



Cf = Ci n E d 



^ 1 r 



Ci n n s^, C out — c ou t n s d 



as shown in figure 1(c) Now, by the Schwartz-reflexivity of vr (see [2J, p. 1159) and /it, 

v(s){vr(s) - /) 



1 
2~7ri 



and therefore 



We write this as 



where 



R(0) -I = U 
R(0) -1 = RM 

rW = n 



-ds 

1 

ni 



1 

2ni 



v(s)(v R (s) 



-ds 



KS)(VR - /) 



s 



(67) 

(68) 

(69) 
(70) 



1 



1 



RW = 



M s ■"fl s - 

90, s 



i? (3) = 3? 



M(a)M*)-i) 



(71) 



10 



Hence using (p3)) . 



If ds 

V>(s){v R (s) - I) — 



TTl .ln+i ir7+ ut 



-f (M S ) -/)-«srw- 

7Ti 7 9G ga s 



[2t-Mt 1/3 -l] 

E lo s(^-i) 

[2i-MW 3 -l] , 1 

E 1 - + lo S7 + 1) + ^ log7 - log (iZa 2 (0) - VT^iMO)) [ 
E 1 -?(--+log(^ 7 J+l) + ^log(- 



1 ' + 2 1 ° g V0 



3.1 Calculation of = 5? 

First, we compute vj^ explicitly. 

Lemma 3.1. For z € 90£, £/ie jump matrix vr(z) can be written as vr = I + -^v q R + Vr where 



1 / rfi/3 2 + ci/T 2 i(di/3 2 - ci/3" 2 ) 

2 Udi/3 2 - ci/T 2 ) -(di/3 2 + Cl /T 2 ) 



5 a 7 

Ci = , dl = 

72' 72 



and 



Proof. On <90£, vr = m p -m p ^_. On this contour to p _ = m*- 5 ' 00 -* and m p + is given by (|60|) . Thus for 



m (6,oo) e -,a«a.^-l [ ( _ ga ( z ) 



2/3 N 



3 r ^ 

2«( z ) 



-i7r/6 1 / 1 — 1 

e /D g e _ _ 



For < arg(s) < 4^, consider (see ([51])) 



v 7 \Ai (s) ^ Ai (w 2 s) y 



with w = e 2i7r / 3 . From Abramowitz and Stegun [T] (10.4.59) and (10.4.61), for | arg(s)| < 7T, 

--(2/3)s 3/2 



Ai( s ) 



1_ * +0 ,' 



wherein ci = ^ and di = — Also note the identity 



72 ■ 



Ai(s) + wAi(ws) + w 2 Ai(w 2 s) = 



11 



and, from Abramowitz and Stegun (10.4.11.13), 

W[Ai(s),M(ws)] = 

W[Ai{s),Ai(u; 2 s)} : 
W[Ai(ws),Ai(w 2 s)] 



-27T/6 



2tt 
1 

1 



s i77/2 



Using det*(s) = W[Ai(s), Ai(w 2 s)] = ^e"^ 6 , we have 

tf _1 (a) = 2tt 

Using |A(z) 3 / 2 = a(z), equations (J77J) and (JTHJ) yield 

Ai(q 2 / 3 A(z)) = 



w 2 Ai'(w 2 s) -Ai(w 2 .s) 



-got 



Ai'( g 2 / 3 A(z)) = 
Ai(<A 2 / 3 A(z)) = 



20F(q 2 /3A)l/4 

(q^A) 1 /^"^ 



2V^ 

J-rr/6 qa 



1-^ + 



1- — +0 

get 



1 



|ga| 2 
/ 1 



20F(g 2 / 3 A) 1 /4 



go 



VlH 2 



^ 2 AiW/ 3 A(z)) = -- 



o; 2 ( g 2 / 3 A) W« / + 4 Q / 1 



(80) 

(81) 
(82) 

(83) 

(84) 

(85) 
(86) 
(87) 



We insert the asymptotics ([84]) - ([87]) into ([83]) resulting in the asymptotic formulas for < arg(g 2 / 3 A(,z)) < tt 
and qa large, 



*- 1 (g 2 / 3 A(z)) 



1 -1 

—i —i 



1 

qa 



d\ -ci 



O 



\qa\ A 



-qa 



ct 3 /6 



It is straightforward to compute an analogous expansion for ^E' _1 (s) for the other values of arg(<7 2 / 3 A(z)) 
in equation (|76|) . To do this one must use the asymptotic formulas ([77]) and ([78]) as well as the additional 
expansions (10.4.60) and (10.4.62) from Abramowitz and Stegun [T]. Namely, for | arg(s)| < 27r/3, 



Ai(-a) 



-1/4 



s l/4 

M(-8) = —= 
\ 7T 



sm 



cos 



2 3/2 7T 

3 4 



3ci 

2s 3 / 2 



COS 



2 3/2 7T 

3 4 



^. s 3/2 4. l\ _ Ml 

2s 3 / 2 



sm 



(89) 



(90) 



After carrying out this computation, the first two terms in the expansion are the same in all four regions. 
In other words, (|88[) is valid not only for < arg(g 2 ' 3 A(z)) < 27r/3 but for all regions in the definition of $ 
in ([75]) . Inserting the expansion in ([55]) . equation ([75]) reduces to 



_ T ,J_f ^p 2 + Cl /r 2 l (d 1 /3 2 - Cl /r 2 ) 

R + 2qa {tidiP 2 - d/T 2 ) -(d^ 2 + Cl p~ 2 ) 



o 



1 



\qa\' 



for all zedO ( . 

Now we explicitly evaluate 



(91) 
□ 
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Lemma 3.2. We have 



89(7-1) 24(J7 
(7-1) 



(7-1) 



1/2' 



1/2 



8g(7-l) 1/2 24g7 



v 8q(7-l) 1 / 2 24 9 7 8g(7-l) " r 24<j7 

Proof. From Lemma 13. 11 it is sufficient to compute the integrals 



Jl 2m 



^ . ds and ^2 = 



1 



27r« iao,; /3(s) 2 a(s)s 



<is. 



(92) 



(93) 



1/2 



We will use the relations £ = e* 6 * and 

• 2 (7-i) 1/2 r/n 7-2 . /0 C \ 1 /0 C \ / 7 -l 

sm( c) = _ , cos( ^ c) = ^ ; sm ^_j = _ cos ^_j = (^__ 

Note that a(z) = |(z — £) 3 / 2 G(z) for an analytic function G(z) in 0£ (see the bottom line at p. 1157 of [2]). 
Hence by residue calculations, 



(94) 



h = 



-dz = 



and 



I? = 



2m J dQ( 2{z - Z)(z - 1/2 G(z)z 2 (£ - £) 1/2 G(0£ 
3(z-0^ 



(95) 



2m J QOt 2{z - 2 G(z)z 



dz 



3 / d 


\{z-l) 1/2 - 




2 [dz 


[ G{z)z _ 


) 



1 _ (C-0 1/2 G'(£) _ (a-s;) 1 / 2 

2(t-o i/2 G(m G^n cm 2 



(96) 



But since a(z) = |(z — £) 3 / 2 G(z) and a'(z) = — ^ ^^~\J ( z — — £), a straightforward computation yields 
that 

a'(z) 7^+1 



G(£) = lim 



-(^?) 1/2 . 



and 



Using (|94|) . we obtain 



and 



Therefore, 



G , (0 = 37Ae + 2 )( e-ov 2 



20 



/1 



2e 2 (e-o i/2 



p(i) 



ID 
22 



and 



^12 — ^21^ — ^ 



4(7 


-1) 


3 




-5T 


(4(7 






72 Jl ~ 





3(7-1; 

57 



1 



1 



" 89(7 - 1) 24g 7 
1 (7-l) 1/2 



8(7(7 - 1) 1/2 24 97 



(97) 
(98) 

(99) 

(100) 

(101) 

(102) 
□ 
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3.2 Bound on = R 



ds 
s 



We begin by establishing the leading term of a(z) for z near £. 

Lemma 3.3. For 1 < q < 2t and for z such that \z — £| < min{ i, |£ — £|}, 



a(*) - - ( 7 - l) 3/4 (^ - e) 3 / 2 e -W 4 e- 3 ^/ 2 



< 



50( 7 - l) 3/4 |z-£| 5/2 



(103) 



Proof. Write z = £(1 + e). Then |e| = \z — £| < min{|, |£ — £|}. Under the change of variables s = £(1 + eu), 
equation (|54p for a(z) becomes 



a(z(e)) 



7 e 



4V£ 



M 1 



£e 



1/2 1 + 



£-£ / + 



Using £ = e l6,<: and (|M|) . we have 



7C 3 / 2 (1 + £)y/£ £ _ _ ^3/4^3^/^3/2 _ / _ ^3/4, _ ^3/2 e -3i7r/4g-3ie c /2_ 



(104) 



(105) 



For the integrand in (fTUll), using the inequalities |(1 + w) 1 / 2 — 1| < |w| for \w\ <land|(l + w) 2 — 1| < 10|w| 
for |io| < |, and using the fact that jY^iy ^ |£jf an< ^ ^ — \£ 2 ~£\ > we °btain 



1 + -^=i 



1/2 i i J^ u 



£-£ / (l + e«) 2 



1 



< 



50|e| 50|z-£| 



ie-£l l£-£l 



Therefore, we obtain (| 10311 . 

Lemma 3.4. For L + 1 < q < 2t - Mt 1 / 3 - 1, there is a constant c > such that 

c{2t) 2 



|i? (2) | < 



q 3 / 2 (2t-q) 5 / 2 ' 



Proof. On <90£, |z — £| = <5|£ — £| < ^j|£ — £|. Hence from Lemma T3.31 we have 

Hz)\ > I( 7 - l) 3 / 4 |z - £| 3 / 2 = ^(7 - l) 3/4 k - CI 372 = ^ 3/2 (2^1 
for 2 G <90£. Therefore, as and ^ are bounded on dO^, 



3/2 



|i? (2) | <c' / 
Jd 



r 



ao t ? 2 (7-l) 3 



ldsl = 



2ttc'7 3 ^I£ - £1 c 7 2 
o 2 ( 7 - l) 3 ~ 2( 7 _ 1)5/2 



for some constants c', c > 0, as |£ — £| 



4(7~1) 1/2 



(106) 
□ 

(107) 
(108) 

(109) 
□ 



3.3 Bound on = R 



Since fi(s) and 1/s are bounded on \R^\ < c'\\v R - I\\ L i (c + } for some constant d > 0. But on C- 7 , 
v R (z) -I = 0(e- 2 i a ^). Hence 

\R^\<c\\e-^\\ LHCt) , (110) 
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for some constant c > 0. For z = e £ C+ (hence 9 C < 9 < ir), using 



we have 



^(e^-O^-O = l(e 40 - e l ^)(e # - e - < *«)| 1/a e <(,r+ * )/a , ( ln ) 

a ( e tf)=_Z/" 1 - e' 9 c)(e^ - e~ i9 °) ■ ie^tty 

, 61 (112) 

f (f\ . 1/2 f ^+ftA . 1/2 ' 

= 7 / cos — sin ' — - — sm ' 



2) \ 2 J \ 2 

(Recall that 7 = y .) Note that a(£) = 0, a(s) is real and positive on C± , and a(e l6 ') increases as increases. 
Lemma 3.5. For 1 < q < 2t, 

a{e ie )>r^VlW^){0-e c f (113) 

/or 6 C < 9 < 7T. 

Proof. We consider two cases separately: # c < § an d # c > ^. 

Start with the case when # c < -| . We consider two sub-cases: 9 < and 9 > =^-. When 9 < 
< § < §, < ^±£<= < I and < < §. Hence using the basic inequalities cos(a;) > \ for < x < § 

and sin(x) > —a; for < x < §, we find that 



a{e w ) > 2. / (0 + fl c )V2(0 _ e c )V2 # > _L / _ W = _L(0 _ 0c f. (114) 



^)>^ _/ (0 + ^) 1/2 (0-0 c ) 1 / 2 #>^ 
When 9 > from the monotonicity of a(e ) and using (|114[) 



a( e ^)>a(e^)>^(^-0 c ) . (115) 



For < 9 C < § and ?f < 9 < n, we have ^ - C > § > §(0 - 9 C ). Therefore, 



a(e«)>^(6-9 c y. (116) 



For the second case, when 9 C > ?, using the change of variables ^ 1 — > 7r 



o( e», =,/;_;• Si „ (|) sin-/* (*±^) (<^i) # . (117) 

Note that < § < f , < 0+(lT 2 ~ 9c) < tt - 6» c < ^ , and < < |. Using the basic inequalities 

sin(x) > for < x < f and sin(x) > for < x < ?f, we find that 

a(e i9 ) > -^D- f ° 0(0 + (vr - # c )) 1/2 ((^ - fl c ) - 0) 1/2 d0 



8^' 



7T — 6 



(118) 

' (TT - C + 2(tt - 0)) (0 - C ) 2 > -^-r(7T - C )(0 - C ) 2 . 



16\/2 V c ; 16-v^tt 

Since = cos (%) = sin (^) < we have 

<*(e i9 ) > ^-jVtCt 3 !)^ - c f. (119) 
Combining HUD , HUH), and (fill))) completes the proof. □ 
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Lemma 3.6. For L + 1 < q < 2t — Alt 1 / 3 — 1, there is a constant c > such that 



i fl,3, i <- ^gf-V (l20) 



Proof. Let e ie * be the endpoint of Cy 4 " on dO^. Note that since radius of dO^ is <$|£ — £| = 45- 



,/W - 1 

6) c > 4j¥J . (121) 



Using Lemma 13751 and changing variables, 



where 

Using the inequality J°° e e£x < -\ for a > 0, 

" " il(Cl+) W7(7-1)V ^" 128<53 9 3/2 (2t _ 9)5 /2- ^ 
Hence from (fTT0|) we obtain ([120]) . □ 

3.4 Bound on i^ 4 > = 3ft 



As before, since on UC^ the functions /x(s) and i are uniformly bounded and Vr{z) — I = 0(e _2l?Q ( 2 )), 

|i? (4) | < c'\\v R - /|| £ i (c£uc + t) < ^H^c+uc+j (125) 

for some constants c', c > 0. 

Lemma 3.7. For L + 1 < q < 2t - Alt 1 / 3 - 1, there is a constant c > such that 

Proof. Let 70 = csc 2 ^). Let <5 4 > be a small positive number defined on p. 1152 of [2]. We estimate 
e~ 2qa in the following three cases separately: (i) 2i(l + & t )~ 1 < q < 2t-AIt 1 / 3 + l, (ii) < q < 2i(l + <S 4 )- 1 
and (hi) L + 1< q < — . 

(i) For 2f(l + S^- 1 <q<2t- Alt 1 / 3 - 1, from (6.37) of 0, there are constants ci, c 2 , c 3 , c 4 , c 5 > 
such that 

/■oo 

H^ILmc+uc+O ^ c i / e-^^dx + c 4 e- c ^. (127) 



C2V7-1 



Using the change of variables y = c^qx 3 , 



< e ^4c 3q (y-l)^ + C4e -c 5q ( 128 ) 



< 



C\ c 4 



3c|clg 2 ( 7 - 1) 5 /2 c 2 ? 2- 
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Since 7 > 1 and 7 — 1 < 1 + 64, we find that 

lle 2 m + + < C7 ' - C[2t? (129) 

for a constant c > 0. 

(ii) For =^ < q < 2i(l + 64)~ 1 , note that the radius of is of O(l). Hence a standard calculation in 

Ricmann-Hilbert steepest-descent analysis shows that for z e ft C^t, 9ff(a(.z)) < — c' for some constant 
d > 0. Since the length of ^(C^ U C+ ut ) is bounded, 

Hp 2 ?"!! j 4. < r-'V- 2 ^" < c " < c "' 72 - C "'( 2 Q 2 ( ,o n) 
H e Hi 1 (c+uc+ t ) S S (2c , )V S ?2(7 _ 1)5/2 - q3/2{2t _ q)5/2 UdUj 

for some constants c" , c'" > 0. 

(iii) Consider the case when L + I < q < —. Then < 6 C < j^. In this case, we make a specific 

choice of and C^ ut : 

C+ = + psm^M-> : 25 < p < 1 j 

C+ t = U + psmceW'-to :25<p< 



a 
1 

sin(6> c - ±7r) 



(131) 



The contours are straight line segments from £ to a point on the positive real axis. (Recall that has 
the radius 5|£ — £| = 25sm9 c .) Now we estimate "R(a(z)) for z e U C^ ut . For z G C^, take the contour 
in (j54|) to be the straight line from £ to z. Then one can check from the geometry that 

< arg(l + s) < 9 C , < arg(s) < 6 C , arg(s - f) = 6 C + -J 
-< arg(s -£)< — - 6» c , arg(ds) = C + — . 

ZD 6 

Therefore the argument of the integrand in §M§ is in [2ir, 2tt + ^ + 26 c ] C [2-7T, 27r + |] since < 8 C < j^. Thus, 
the cosine of the argument is greater than or equal to cos( : |) = \. Therefore, for z = £ + psmO c e l ^ 9c+ ^^ G 
C^T, by the change of variables s = £ + y sm9 c e' l< - e ' :+ ^ 7T \ 



Vt(a(z)) < -J 



l + s 



c2 



^(*-f)(*-0 



9 ^ 11 1 yBine e i(^+m> (133) 

7Bin 2 ^cos(^) /* J 1 + 7^T) e '\ .I i{ec+W j2 



l\l 2 Jo 1 1 + y sin y c e 6 I 

Using the inequality |1 + xe^| > | sin</>| for all a; e R, and using < 9 C < and |y| < _ sin( -g + 2^ ) 5; 2, we 
have 

7sin 2 6» c cos(^) f p 



^lzT)"> (131) 



27 V 7 



17 



for z € C^t. For z e C£ ut , taking the contour in (fi>4")) to be the straight line from £ to z, we can check that 
< arg(l + s) < 6 C , < arg(s) < 6 C , arg(s - £) = i 

arg(ds) = 6» c - -. 



7T 

6 



| - C < arg(s - £) < |, 



(135) 



Hence the argument of the integrand in (|54|) is in [— 6 C — ^,29 c ] C [— §,§]• Therefore, for z = £ 
psin^e 1 ^"^) g C+ 



U(a(z)) < 



< 



< 



7 
4\/2 

7 sin 



1 + s 



(*-0(*-0 



\ds\ 



|1 + yshi9 c e 



-i%\2 



1 — i— ye 



^c-f)|V2 



1 /7-l\ 



3/2 



(136) 



< -- 



27V2V 7 / 
1 f7-l/j 



n 3/2 



(2<S) 3 / 2 . 

27v / 2\ 1 J 
From (|134[) and (|136[) . arguing as in (|130[) . we obtain 



3/2 



P 2ga| 



< 



c"(2t) 2 



for a constant c" > 0. Hence we obtain the estimate for |i?( 4 )|. 



(137) 
□ 



3.5 Bound on = K 



1 qotds 



-V 



qa 



R 



Lemma 3.8. For L + 1 < q < 2t - Mt 1 ^ - 1. i/iere is a constant c > smc/i t/iai 



|i? (5) | < 



c(2i) 2 



g 3 / 2 (2i-g) 5 /2- 

Proof. As — / = (1 — £r) -1 Cr/, and as (1 — C) _1 and £_ are uniformly bounded. 

\\H - I\\l*<p) < cq\\£rI\\l*(t>) = co\\£-(vr - < ci\\v R - 



(138) 



(139) 



for some constants Co, c\ > when t is large enough. Below, we assume that t is large enough so that the 
above estimate holds. Now 



|i? (5) | < 



qa 



<2|| At -7|| i2(a ^ ) 



\ds\ 



< 



2ci\\v R - I\\ L '(e) 



qa 

qix 



R 



qa 



L 2 (dO ( ) 



L2(ao 5 ) 



(140) 



Since /3(z) = (■f^) 1 ^ 4 is bounded above and below for z € 0£, v q R a (z) in Lemma [3.11 is bounded. 



Using (|108| and the fact that the radius of is 6\£ — £| 



7, _ iS^P^T 



we have 



qa 



< 



C27 



c 2 (2t) 2 



9 2 ( 7 -l) 5 /2 g 3/2(2f-g)5/2 



(141) 
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for a constant C2 > 0. 

Write \\vr- I\\ 2 L 2 { ^ } = \\v R - I\\ 2 L2{ao) + \\v R - I\\ 2 L 2 { ^ aoy As v R - I = ^- + in Lemma O is 
bounded by t^j- for a constant c 3 > 0, we have, as in (|14ip . 



UvR-iWhw = 2|K-/|li 2(90)4 < g3 / 2 ( 4 2 f! ) g)5/2 ( 142 ) 

for a constant C4 > 0. On the other hand, 

IK - I\\h(E\ao) = 2 IK ~ 7 H^(c+) + 2 IK _ / "i 2 (c+uc+ t ) 

< c3(||e- 29Q || L2(Ci+) + ||e 2 ^|| i2(CituCo+ut) ) (143) 

I 1 I I _2qrcK I I \ 



<C3(||e- 2 ««|| /jlf ^ + ||e 2 ^|| J - irf7+uf7+ , 



for a constant c 3 > 0, since e~ 2qa < 1 for 2 6 and e 2?Q < 1 for z G C*+ U C+ Ut . Hence from (fl24j) 
and (|130p . we have 

ll^-/||| 2(S \aa)< g3 / 2 ( 4 2 f! )2 g)5 / 2 ( 144 ) 

for a constant c 4 > 0. By combining ([Ml]) . fliS]) . and (|T24"|) . we obtain ([T38]) . □ 

3.6 The Airy part 

From Lemmas E3 EH E3 and[3?8l we find that, for L + 1 < q < [2t - Mt 1 / 3 - 1], there is a constant c> 
such that 

„ c(2t) 2 c(2t) 2 



and 



i= 2 

We need the following result. 
Lemma 3.9. We have 

pM-Mt 1 '"-!] 2 

lim limsup V ^ J =0 (146) 

L+l v 7 

[2*-M* 1/3 "l] 2 

lim limsup £ { J = 0. (147) 

q — -L+l 

Proof. We use the following basic inequality. Let a, b be integers. Let s(.t) be a positive differentiable 
function in an interval [a — 1, b + 1] and there is c £ (a, 6) such that s'(x) < for x G [a — 1, c) and s'(x) > 
for a; G (c, 6 + 1]. Then 

X] s(g) < / s(a;)da;. (148) 

q=a "' a - 1 

As a function of < q < 2t, q 3/2^ 2 t_ q ) 5 / 2 decreases for < q < and then increases for |£ < q < 2t. 
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Hence 



[2i-A/t 1/3 -l] 



(2t) 2 



,.2i-Mt 1/3 



< lim lim sup 



q 3 / 2 {2t- q) 5 / 2 



dq 



(149) 



lim lim sup 



4(-3t 2 + 6i<?-2g 2 ) 
3i(2i - qfl 2 q x l 2 



2t-Mt 1/3 



As a function of < g < 2i, 



(2t) 2 



9 2 (2t-g) ; 



decreases for < g < t and then increases for t < q < 2t. Hence 



[2t-MiV3_l] 



< lim lim 

L — >oo t — >oo 



lim lim 

L— >oo t— >oo 



2t-Mt 1/3 



{2tf 



q 2 (2t-q) 2 



dq 



-2t + 2q 1 (2t-q 
q(2t -q)~t ° g 



2t-Mt 1/3 



= 0. 



Now we prove the main result of Section El 
Lemma 3.10 (Airy part). We have 



lim lim 

L,M—>oc t— >oo 



[2t-Mt 1/3 - 


11 / 


E 


log(*«-l) - it 2 


sr=L+i 








--AI 3 


log M H lo 


12 


8 S 24 



12' 



Proof. We first prove that 

[2i-A/t 1/3 -l] 



[2t-Mt 1/3 -l] 



lim lim 

£,M->oo t 



f / 2i 2t 2t 

im ^ Iog^)- E -<K--+log(-)+l) + -log(-) 

^ q =L+i q=L+i \ y y y 



l 



(150) 



□ 



(151) 



(152) 



1 

+ 8(2i - q) + Y2q 

Using 1 + P$(0) - y/j^~TR { 2 \ ] (0) = 1 - 8(2t x _ g) - > | and using Lemmas EH EH EI3 and EE 
we find that 



£-= 2 (^(o) -V^t4i(o)) 



i + ^(o)-Vt^^(o) 



i 

< - 

- 2 



(153) 
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for L + 1 < q < [2t - Alt 1 / 3 - 1], when we take t large enough. Hence using | log(l + x)| < 2\x\ for |x| > h, 



hm hm 

L,M— >oo t^oo 

q=L+l 

[2t-Mt 1/3 -l] 

< lim lim 4 > 

q=L + l 
[2t-Mt 1/3 -l] 



i + ^^-Vt^WAo) 



i=2 



, (2t) 2 (2i) 2 
< lim hm 4c > — 1 — 

~ L.M^oot^oo q 3 /2(2t-q) 5 /* q 2 (2t-qy 



(154) 



using (|145p and Lemma [3?9l 

On the other hand, from Lemma 



[2t-Mt 1/3 -l] 

£ io g (i + i?( 1 2 ) (o)- v ^T^ 1 1 ) (o)) 

Using -x 2 < log(l - x) +x < for < x < h 



[2i-A/t 1/3 -l] 

E lo s ( 1 



1 



8(2t-q) 12q 



[2t-Mt 1/3 -l] 

E M 1 -^: 

q=L + l 
[2t-Mt 1/3 -l] 



8(2t-q) 12q 



. [2t-Mt 1/3 -l] 

T^)+ E 



8(2* - g) 12? 



< E 

q=L+l 
r2t-Mt 1/S 

< 



+ 



64(2t-q) 2 48(2t-q)q UAq 2 



1 1, 



64(2i-g) 2 96t K 2t-q q' 144g 5 







(156) 



as t — ► oo by evaluating the integral explicitly. 

Using (fl"54)) and (fl"56)) in J72J), we obtain (fl"52|) . 

Now we compute each term of the sum in (fl52]l . Note that [2t - Aft 1 / 3 - 1] = 2t - Alt 1 / 3 - 1 - e 
for < e < 1. We have 



[2t-Mt 1/3 -l] 

E 

and 



2< 

-« + 1 

q 



[2t-Mt 1/3 -l] 

^ (2t-q) = -(2<-L-l)(2i-L)--(Aft 1/3 +e)(A« 1/3 + l + e) 

(157) 



[2t-Mt 1/3 -l] 

£ (-g + -) log(2i) = --((2t - Alt 1 / 3 -l-ef- L 2 ) log(2t) 



(158) 



Since for positive integer m 



fclogfc = mlog((m - 1)!) - V] log(fc!) = mlog(r(m)) - log G(m + 1), 



(159) 



fe=i 



fc=i 
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we find that 

[2t~Mt 1/3 -l] 



E 

q=L+l 



- ) logg = ( 2t - Mt 1/3 - e - ^ ] logT(2t - Mt 1/3 - e) - logG(2i - Mt 1/3 - e+ 1) 



L+- )lo g r(L + l) + logG(L + 2) 



(160) 



Note that from Stirling's formula for the Gamma function and the asymptotes (|48j) for the Barnes G-function, 

as z — > oo, 

z + M log r(z + 1) - log G(z + 2) = (L. - £j l 0gz - Z - + £ _ I log(27r) + -L - C '(-l) + o(l). (161) 



Now using the fact that 



K 1 

lim I log A' — 7I = 0, 



AT— +OO 1 (7 

9=1 



where 7 is Euler's constant, we obtain 

[2t-Mt 1/3 -l] 

E 



lim 

t — >oo 



q=L+l 



whj-l logi2t)+1 s log{Mtl/3) 



= 



and 



lim 

t — >oo 



[2*-A/t 1/3 -l] 

E l2^-l2 l0g(2<) -12 7 ' 12* 



1 L 1 

E 



Therefore, using {TS7J), (JTSHP, (flBO]) . |gB5[| , fTM]) . and (fT6T]) . we obtain 



lim 

£ — >oo 



1 



1 



[2t-Af{ 1/3 -l] , 

£ (-'(-7 + ^(7) + 1 ) + 5 be <7 ) + 8(2 t -,) ' 12, 



t 2 - 2tL + — \og(2t) 
_1m 3 - 1 logM + -Llog2 - i log(27r) + 1 - C'(-l) 

i(L + 1) (l, + lo g r(L + 1) + log G(L + 2) + ^7 _ q 



L 1 



Hence using (|161j) and (|162[) again, we obtain (|15ip . 



(162) 



(163) 



(164) 



(165) 



□ 



4 Proof of ([9]) in Theorem ll.lt computation of F(x) 

Recall equation (|37[) which we rewrite here: 



[2t-Mt 1/3 -l] [2t+xt 1/3 ] 

log(e- i2 £> n ) = -i 2 + log(i? i )+ ^ log(A t -_ 2 1 )+ ]T log^) 

g=-L+l g=[2t-MiV3] 



exact part 



Airy part 



Painlcvc part 



(166) 
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For the Painleve part, from [2], 

[2t+xt 1/3 

lim 

t 



/X 
R(y)dy 
-M 



q^fct-Mt 1 



1 o 1 

— M + -logM. 
12 8 8 



By combining (fW|) . (fT5"Tj) . and we obtain 



lim lim log(e~* D„(t)) 

M.L— >oo t^oo 



Let 



i — >oo 



(167) 



(168) 



(169) 



= - £ 7 - -/ + Qdy + - 1 log |x| + 1 log 2 + C'(-l). 

5 Proof of (1101) in Theorem computation of E(x) 

Ij{2t) = — f e 2tcos8 e ij6 d6. 

27r J -It 

Set (see g]) 

L>++(t) = det(Vfc(2t) + /i+fc+a(2*))o<j,fc<<-i, (170) 
I?7 + (t) = det(/,_ fe (2i) + / J -+*+i(2*))o<j,fc</-i. (171) 
It is shown in Corollary 7.2 of [5] that for 

i=[t + f lf % (172) 
where a; lies in a compact subset of K, we have 

lim \e- t2/2 D++(t) - F(x)E(x)\ = 0, (173) 

t — >CO 

lim | e - t2/2 -*D7 + (t) - F(x)E(s)| = 0. (174) 



In [5], the above results are shown (with x in place of x) for the alternate scaling t = £ — ^i 1 ^ 3 . With the 
above scaling £ = [t + ft 1 / 3 ], we find x = x(l + ft" 2 / 3 ) -1 / 3 . Since x is in a compact set, and E(x) and F(x) 
are continuous, the above results follow. From equations (|173p and (] 1 T4[) . for a fixed x g K, 

F(x) 2 £(a-) 2 = lim e- t2 - t £>t+£>7 + , £ = [t + §i 1/3 ]. (175) 

Let 7Tj(z]t) be the monic orthogonal polynomial of degree k with respect to the weight ^e 2tcos d9 
on the unit circle, as introduced in Section [TJ It is shown in Corollary 2.7 of [3] that (cf. (|5zT)) above) 

oo oo 2 

e -* 2 / 2 i?+ + (t) = n 4 +2 ( t )(i - 7T 2J+2 (o ; *)) = n , , K2j " +i , n , v ( 176 ) 

^ ~5 l + 7T2i+2(0, t) 

oo oo 2 

e -t V 2-t D - +it) = Y[4 { t)(l + n 2j+1 (Q;t)) = J] z "^nTT' (177) 
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where the last equalities in (| 176[) and (|177|) use the basic identity (see e.g |23j ) 

l-^(0;t) = 



K k-1 



Using equations (|176[) and (|177p , we can write 

= II 7F+ = ^-i II ^ 2 -i(*)(! + ^(0;t)), 

I D -+ t 

D^ + = D- L + H -J^=D- L + H K^ a _ a (t)(l-7riy_i(0;t)), 

and hence we have 

2£-2 £-1 

D+\DJ + = D+\Dl+ [J K - 2 (0.n[(l + 7r 2j (0;0)(l-7r 2j+1 (0;0)]. 

fc=2L-l j=L 

Using (recall that D( = dct(Ij-k(2t))o<j,k<e-i is the £ x £ Toeplitz determinant given in j27 



a-l 

3=0 



we find that 



= £>++ U-+^zl [] [(1 + ^.(0; t))(l - 7r 2j+1 (0; t))]. 

Inserting this equation into (|175p . and using (|29l) . 

£>++ £T+ c_1 

F(x) 2 E(x) 2 = F 2 (x) lim e~« ^ J] K 1 + ^(0;t))(l - 7^+i(0;f))]. 
Hence we find (cf. (|36|) ) 

£(x) 2 = lim e -t D ^-i^ + J][(i +7ra .(0;t))(l - 7r 2j+1 (0;t))], 1= \t+U^ 

In analogy to equation (|37[) . we write 

/ ^++^-+ [ *-^ tl/31 

21og£(x) = lim lim - t + log L ~ 1 £ + V log[(l + tt 2j (0; t))(l - ^- +1 (0; t))] 

— *-oo t— »oo \ J-^2L — 1 — 



Exact part 



Airy part 

[t+t* 1/3 -i] 

£ log[(l + 7r 2j (0;i))(l-7r 2j+ i(0;i))] 
J=[*-f * 1/3 +i] 



Painlcvc part 



We compute each term as in the case of log-F(x). 
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Lemma 5.1. (The Painlcve part) We have for x < and M > 0, 



[ t+ | i V3_ 1] 

Mm log[(l + 7r 2j (0;t))(l-7r 2j+1 (0;i))] 

""°°i=[t-¥ * 1/3 +i] 



9(y) 



— M 



v^2 M 3/ 2 _ v / 2(~^) 3/2 



(187) 



Proof. This follows from the following results in Corollaries 7.1 of [5]: for x in a compact subset of R, 



lim 

£ — >oo 



^ log(l + 7r 2j - +2 (0;t)) + \ogE{x) 

j=[i+ftV3] 



lim 

£ — »oo 



= 0, 



0. 



5^ log(l-7r 2j+1 (0;i))+log£(; 

3=[i+ftV3] 

We remark that in [5], the notations i>(x) = — R(x) and u(a;) = — q(x) are used 
Lemma 5.2. (The Airy part) 



lim lim 

L,M— >oo t^oo 



[*-^ 1/3 ] 



£ log[(l + ^(0; t))(l - 7r 2j+1 (0; i))] - f t - y^M^ - ( 2L - i) log 2 



(188) 

(189) 
□ 

0. (190) 



Proof. Using ([64 

[t-M t i/3] 



2[t-f-t 1/3 ] + l 



log[(l+7r 2j (0))(l-7r 2i+1 (0))]= £ io g ('l + ./^3i2 11 (0; g )-^ J R 12 (0;g)y (191) 

j=L q=2L V V 7 V7 / 



Using (|70|) and the fact that ^7 — LR^(0) — R\^(0) = 0, which follows from Lemma I3~2l this equals 



2[t-M t 1 /3 ] + 1 
q=2L 



7-1 



2 [t_M t l/3] + 1 



V log ( 1 + St 2 (V7^i4? (0; <?) - i$ (0; ?)) 



q=2L 



V7 + Vt^t 



(192) 



From LcmmasEHEHO and^Ei the same estimate as in (|T35]) holds for Ei^CVT^Mi (°5 <?)-#i2 (°5 ?)) 
for L + 1 < q < [2t - Mt 1 / 3 - 1]. Therefore the same argument as in (|154p implies that the second sum 
in (|192p vanishes in the limit. Therefore, 



lim lim 

L,M— >oo i— >oo 



[i-M* 1 / 3 ] 2[i-ft 1 ' 3 ]+l 

2 log[(l+7r 2j -(0;t))(l-7r 2j - +1 (0;t))]- £ log ( 1 

j=L q=2L 



2t-q 



2t 



0. (193) 



We use the Euler-Maclaurin summation formula 



E /(*) = / f( x ) dx + 



f(a)+f(b) f'(b)-b'{a) 
2 6-4! 



+ Err, |Err| < J \f<®(x)\dx. (194) 



Set f{q) = log ( 1 + V ) and write [f - f-t 1 / 3 ] = t - ^-t 1 / 3 - e where < e < 1. Then 

lim /(2*-MtV3- 2e+1)+/(2L ) = 1 bg2; 
t— >oo 2 2 



(195) 
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and 



f(2t-Mt 1 / 3 -2e+l) + f(2L) 

hm — — = 0. 

f^oo 6 • 4! 



(196) 



Also using / (3) (g) > 0, 

2 /.2t-A« 1/3 -2e+l 



|Err| < 



(2tt)' 
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l/ (3) (?)l^ = (/" (2* - Mi 1 ' 3 _ 2e + 1) - /" (2L)) - 0, (197) 



as i ^ oo. Finally, changing variables to w = 1 + y and setting oi = 1 + y and S 2 = y MltH^i 1 ; 

r 2i-Mt 1/3 -2e+l / 



2L 



2i-g 



log 1 + W \dq = At I (1 — w) log(w)dw 



2t 



l+d 2 



Si 



At 



1 



1 



if — -ur logiu + — ur — w 



(198) 



= i - ^M 3/2 - 2L log 2 + O f-i 



Hence the result follows. 

Lemma 5.3. (The exact part) We have 



lim lim 

L— >oo i— >oo 



log n" 1 L -(2L-l)log2 

^2L-l 



□ 



(199) 



Proof. Note (see [J) that 



= det 



= det 
= det 



— ^ [ e i( j~ k ^ 9 - e i(0+k+2)8^ e 2tcose ( j i Q 



0<j,k<L-2 



— f e l{:i+1)8 sm{{k + l)8)e 2tcose dd 



(200) 



2 sin((j + 1)0) S m((k + 1)6), 2t 



0<j,k<L-2 

sin 2 9e 2tcos9 d6 



7T Ir, sine 



sin ( 



Changing variables to x = cos9 and noting that S£l£LhMff) j s a polynomial in a; of the form 2 J x J + • • • (a 
constant multiple of the Chebyshev polynomial of the second kind), we have 

- 2 j+fc+i rl 



0<j,k<L-2 



Dili = det 



/ x J+k ^l-x 2 e 2tx dx 

* J-l 



0<j,k<L-2 



2 \ L ~ 1 2 ( L - 1 )( i -2) 

n) (L-l)! 
A steepest descent analysis yields that 

e 2t(L-l) 



L-l 



n - ^i 2 n v 1 (/ -'v ( 2o1 ) 

l<k<t<L-l j=l 



lim £>++ 



L-l 



^-i )2+( ,-i) /27r( ,-i) (£ _ 1}! Jm[ ., ji^ i» - «i a • n vwe--^- ] = 1. 



(202) 
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The multiple integral is another Selberg integral (corresponding to a Laguerre ensemble) which is evaluated 
explicitly as (see (17.6.5) of [2"2"] 1 



L-l 



G(L + i)G(L) 



\<k<KL-\ 



Hence 



2t(£-l) 



G(§) 

G(L+l)G(L) \' 1 _ i 



lim D + + • 

^(L-l)2 + ( L-i)/ 27r (L-i) " G( 3) J 



(203) 



(204) 



Similarly, 



det 


2ir] 


det 


V. 


det 




det 





f e i{j-k)6 + e i{j+k+\)6 \ g 2t cos e d Q 



0<j,k<L-l 



3 2tcos9 d9 



cos j H — 9 cos fe 



1 c os((j + 1)0) cos((fc + 1)0) 

COS I COS £ 



0<j',fc<£-l 

e 2tcos9 d9 



(205) 



0<i,fc<L-l 



cos^ I - I e 2tcose d9 



0<j,k<L-l 



Setting .x = cos0 and noting that — is a constant multiple of a Jacobi polynomial pj 2 ' 2 \x) of the 
form 2 J x J + ■ ■ ■ , we have 



2 L(L-1) 



det 



1-x 



0<j,k<L—l 



oi(i-l) /■ ^ 

W/^ n i 2 rw 

■'L- 1 - 1 ! l<k<i<L 3=1 V 



1 + X 



(206) 



A steepest-descent analysis implies that 

p ith 

lim L>; 



j-L 2 -L/2 



The multiple integral is also a Selberg integral (see (17.6.5) of [12]), and we obtain 



lim DV+ ■ 



„2iL 



AS. L \t L2 - L / 2 TT L 

Using pWj) . PP5|) . and 



G{L + l)G{L+\) 
G(i) 



= 1. 



lim D 2 l-i 

t — >oo 



„4tL-2« 



X2t) 2L2 - 2L +i(2n) L - 
which follows from Section [2] (the exact part for F 2 (x)), we obtain 



-G(2L)) =1, 



lim 

t — »oo 



log L n X -log 



2 2L " L [G(L+i)] 2 G(L)G(L + l) 



r 1 



G(i)G(|)G(2i) 



The result is now proved using the properties (|46[) and (|4"8")l for the Barnes G-function. 



(207) 

(208) 

(209) 

(210) 
□ 
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Combining equation (|186[) and Lemmas 15.11 15.21 and 15.31 we obtain 

2\ogE(x) = jf^ (q(y) dy ^t^i!l _ I log 2. (211) 
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